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Abstract. We study the magnetic and superconducting instabilities of the 
periodic Anderson model with infinite Coulomb repulsion U in the random phase 
approximation. The Neel temperature and the superconducting critical temperature 
are obtained as functions of electronic density (chemical pressure) and hybridization V 
(pressure) . It is found that close to the region where the system exhibits magnetic order 
the critical temperature T c is much smaller than the Neel temperature, in qualitative 
agreement with some Tn/T c ratios found for some heavy-fermion materials. In our 
study, the magnetic and superconducting physical behaviour of the system has its origin 
in the fluctuating boson fields implementing the infinite on-site Coulomb repulsion 
among the /— electrons. 
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1. Introduction 

The superconducting and magnetic properties of heavy-fermion materials have recently 
attracted much attention because of their non- conventional character. |TJ |2] These 
materials have very large specific heat coefficients 7, indicating very large effective quasi- 
particle masses, hence the designation heavy fermions. Some of these materials, order 
antiferromagnetically at low temperatures (examples are UAgCu^, UCut, XJ^Znyj) while 
others (such as UBeiz, CeCu 2 Si2, UPt s ) order in a superconducting state and others 
show no ordering (such as CeAl 3 , UAuPt^, CeCu & , UAl 2 ). [I] Some compounds exhibit 
phases where antiferromagnetic order coexists with unconventional superconductivity. 
Examples are: UPd 2 Al 3 (T N = 14.3 K and T c = 2 K), UNi 2 Al 3 (T N = 4.5 K and 
T c = 1.2 K), CePd 2 Si 2 (T N ~ 10 K and T c ~ 0.5 K) and Celn 3 (T N ~ 10 K and T c ~ 
0.15 K). In the prototype heavy-fermion system Cea;Cu2Si2 the coexistence of d— wave 
superconductivity and magnetic order was clearly identified in a small range x values 
around x ~ 0.99. [H] 
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Systems that exhibit both superconductivity and antiferromagnetism at low 
temperature have ratios between the Neel temperature T N and the superconducting 
critical temperature T c of the order of Tjy/T c ~ 1 — 100. The coexistence of both types 
of order can be tuned by external parameters such as external pressure or changes in 
the stoichiometry. Ej 

A description of the normal state properties of the heavy-fermion systems has been 
attempted assuming a generalization of the impurity Anderson model to the lattice. 
Ej In the Anderson lattice the energy of a single electron in an /—orbital (e. g. 4/ 1 ) 
is eo, and the energy of two electrons in the same /—orbital (4/ 2 ) is 2eo + U, where U 
is the on-site Coulomb repulsion. The energy of the 4/ 2 state is much larger than the 
energy of the 4/ 1 state. Thus, if the charge fluctuations at the /—orbital are small, the 
4/ 1 electron may behave local moment. 

The complexity of heavy-fermion systems arises from the interplay between Kondo 
screening of local moments, the antiferromagnetic (RKKY) interaction between the 
moments and the superconducting correlations between the heavy quasi-particles. The 
local moments form in partially filled / shells of Ce and U ions. The absence of 
magnetic order in some cases could perharps be due to complete Kondo screening (below 
the Kondo temperature T#) or to a spin liquid arrangement of the local moments. 
In the normal non-magnetic state the Anderson lattice model predicts Fermi-liquid 
like behaviour and explains the main features at low temperatures, such as the large 
effective masses and the Kondo resonance near the Fermi level. But the main technical 
difficulty is the competition between the Kondo compensation of the localized spins 
and the magnetic interaction between them. This interaction is mediated by the 
conduction electrons (RKKY-type). Related to this competition is the effectiveness 
of the compensating cloud around each /-site. The size of this cloud has been subject 
of controversy. While some arguments show that it should be a large scale of the order 
of vf/Tk [U, other arguments claim to be ~ a (a is the lattice constant). JHj This is a 
relevant issue and is related to Nozieres exhaustion problem which states that there are 
not enough conduction electrons to screen the /-moments. 

It has been proposed that the mechanism for superconductivity lies in the strong 
Coulomb interaction between the /-electrons, not in a phonon mediated attraction. 
Using Coleman's [H] slave boson formalism together with a large-iV approach, various 
attempts have been made to search for the existence of an effective interaction which 
might be responsible for superconductivity in the infinite-f/ Anderson-lattice model. 
It was proposed JUj that slave boson fluctuations can provide an effective attraction 
between the electrons to leading order in 1/N. Later, a calculation of the electron- 
electron scattering amplitude to order 1/N 2 revealed an effective attractive interaction in 
the p and d channels, which was interpreted as a manifestation of the RKKY interaction, 
showing that spin fluctuations are an important mechanism. JTI] The inclusion of /°, 
f 1 and f 2 states, using two sets of slave bosons, was also considered in the context of 
the Anderson lattice as a possible description of high-T c superconductors. [T2] 

The magnetic order in the ground state of Kondo Insulators has been studied by 
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Dorin and Schlottmann in the framework of the Anderson-lattice model 13.. The same 
authors have later studied the effect of orbital degeneracy and finite [Zona ferromagnetic 
ground state (their approach did not generate RKKY interactions, thus preventing the 
study of antiferromagnetic order) . [T4*] 

In this work we consider the slave-boson approach to the infinite- U Anderson- 
lattice model. We treat the boson fields at the mean-field level, thereby enforcing the 
constrain of one /-electron (at most) per site only on the average. By splitting the boson 
operator into a condensate part and an above-the-condensate term, which describes 
fluctuations, we compute the magnetic and pairing susceptibilities at the random-phase 
approximation (RPA) level. For spin 1/2 particles the condensate density at moderate 
temperatures does not change much relative to its ground-state value. Therefore we 
do not expect our results to be of less quality then those characterizing the ground- 
state properties. We search for the critical temperatures (T/v and T c ) at which anti- 
ferromagnetic order or superconducting (s or d-wave) pairing occurs in a normal non- 
magnetic system. We find that the value of T c is much smaller than the magnetic 
temperature T N . Unlike T N , the superconducting temperature monotonically increases 
with externally applied pressure. 



2. The model and the RPA solution 



The PAM Hamiltonian is given by 
where 



H = H C +H° f + H cf + Hu, (1) 



H° f = J2(e -ti)flfi,*, (2) 

i,<7 

Hi = £(«% - »H,A > ( 3 ) 

k.a 

H cf = VY.(c\j ha + flc ha ) , (4) 

i,cr 

Hu =U^2 rii^n u . (5) 

i, 

The c and / operators are fermionic and obey the usual anti-commutation relations. 
The hybridization potential V is assumed to be momentum independent. The term Hu 
represents the strong on-site repulsion between the /-orbitals. We consider U = oo. We 
implement the condition U = oo within the slave-boson formulation due to Coleman, [H] 
in which the empty /-site is represented by a slave boson 6j and the physical operator /, 
in equation is replaced with b\fi together with the constrains of only one /—electron 
per site. The implementation of this constraint amounts to introducing a Lagrange 
multiplier A which will renormalize the bare /—level energy from eo to tf = €q + A. We 
split the boson operators in two terms 

&+ = ViVi<W+4, (6) 
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where z represents the boson condensate and Bl represents the fluctuations above the 
condensate. This procedure leads in leading order to a mean field Hamiltonian [TBI I16j . 
The corresponding mean field equations can be written in terms of the Fourier transform 
of the Green's functions 

GffAk, t-t')=- (T T f^(r)fl(r')) , (7) 

GccAK r-r') = - (T T%(T {r)c\y)) , (8) 

GcfAl r - t') = - {T T%a {r)fly)) , (9) 



as 



T 

Y,Y,GffA k > iuj n) , (io) 



TV 

6 fc,cr 



and 



VT 

e f = £ °~ -7z^rY.Y. G cfA k i iuJ n) , (ii) 

V zly s r _ ^„ 



where N s denotes the number of lattice sites. Equation (jl(J|) states that the mean number 
of electrons at an /-site is nf = 1 — z. For a given number of particles per site, n, these 
equations must be supplemented with the particle conservation condition which yields 
the chemical potential fi for any temperature: 

T 

n = 1 - z + —Y,Y. G ccA k > iuJ n) ■ (12) 
The fluctuations beyond the mean field approach are described by the Hamiltonian 

k,q,cr 

and will be considered in the calculation of the magnetic susceptibility and 
superconducting correlation functions below. The calculation, even at the RPA level, 
of the correlation functions requires the knowledge of the boson propagator. The full 
calculation of the latter is a technically difficult problem by itself, and is still unsolved. 
There are, however, 1/N calculations of D( k, r — r'). [EIE31E] Here we follow the work 
of Evans jTTj and use an asymptotic form for the boson propagator given by 

D(K r - t>) = {T T B^{r)B\y)) ~ I . (14) 

We also adopt the same approximation for the propagator D(k,r — r') = 
{T t bI {t)B^ a {r')) . In the calculation below we shall use mean field fermionic 
propagators. 

The transverse spin susceptibility for the / electrons is defined as 

X_+(g> n ) f P dre^(T T S'(q,r)S + (q,0)), (15) 

Jo 
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where (3 = 1/T is the inverse temperature, T T is the chronological order operator (in 
imaginary time), S~(q) = Y^pfpifp+q,^ an d S + (q) = [S~(q)]^. The calculation at the 
RPA level yields 

/ ( ~* ' \ T f / f (q,iu; n )[l - JT C / C (q,iuj n )] 

y\ _{q,Ujj n ) = =-r 1 =~rr = , 16 

[1 " JTfM H] 2 - J*T f / f (q, iu n )V%(q, iu n ) 

where J = V 2 /(NX). The result f)16|) holds for all values of rif and is a generalization 
of that obtained by Evans [T71 [TH] for the case rif — 1. The functions T(q,iu n ) above 
are given by: 

T f / f (q,iu n ) = -— E Gf f (p,iu m )Gff(p + q,iu m + iu n ) , 



p.lLOn 



Tl c c (q,iu n ) = - — E G cc (p,iu m )G cc (p + q,iu} m + iu n ) , 

^ p,iu m 
1 

Tf c (q,iu n ) = - — E G c f(p,iuJ m )G fc (p + q,iu m + iuj n ) 

There are three possible superconducting pairing susceptibilities that one can define. 
These refer to Cooper pairs of either c— electrons or /—electrons, and a hybrid Cooper 
pair with a c— and an /—electron. We consider the correlation function: 

A dd (q,iu n ) = /V- E r ? (fcO^(4)(T T ^ ia (r)d_g i+ . T (r)4 2 i rf t _ A . 2+ . T ) , (17) 

fel,fe2 

where d = c, / and T](k) is the Cooper pair structure factor, assumed to be either 
extended s— wave or d— wave. The hybrid pairing correlation function is defined as 

= [ E (^4 4 (r)c_, i+r!T (r)4 2ji /t ^) . (18) 

This definition has been used previously in a mean field study of the Kondo lattice |19j . 
At the RPA level the Cooper pair correlation function (fTTj) is given by 



1 - J[T J cf (q, iu n ) + Tit (q, iu n )\ 

and the function (|18|) is given by 

A cf {q, iuj n ) = r'J{q, iu n ) + ' , -^t-jyjz . , 12 • (20) 

[1 - JT J fc (q, iu n )\ 2 + J 2 [T J ct (q, iu n )\ 2 

The r(q,iu n ) functions appearing in the previous expressions are given by 
T^(q,iu n ) =— E V 2 (p)G cc (p, iu m )G cc (-p + q, -iu m + iu n ) , 

Tf c (q, iu n ) = - E v(p)G cc (p, iu m )G fc (-p + q, -ioJ m + iuJ n ) , 

r p,iuj m 

T f J(q, iu n ) = 4 E G ff(P> iu m )G cc (-p + q, -iu m + iu n ) , 

r p,iuj m 
1 

T J c f(q, iu n ) = - E G fc {p, iu m )G cf {-p + q, -iu m + iu n ) , 
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3. Superconducting and magnetic instabilities 

The magnetic and superconducting instabilities of the system are signaled by the poles 
of the corresponding susceptibilities. Therefore, we search for the temperature T at 
which the denominators in the RPA expressions for the susceptibilities vanish: 

K m {Q, 0) = [1 - Jf C / C (Q, 0)] 2 - J 2 f f / f (Q, 0)f-(Q, 0) , (21) 

k m (o, o) = i - j[rg(o, o) + r£f (o, o)] , (22) 
K cf (o, o) = [1 - jr£(o, o)] 2 + J 2 [rf/(o, o)] 2 , (23) 

where Q = (n, it, n) and K m (Q, 0), Kdd(0, 0), and K c f(0, 0) are the Stoner factors of the 
correlation functions (|TB|) . (fTU|) . and (J2*U|) . respectively. Since heavy-fermion materials 
are antiferromagnetic materials we seek for poles of K m (Q, 0) at the antiferromagnetic 
wave vector Q. From the definitions of the T(q,iuj n ) functions we see that the 
Cooper pair structure factor r](p) does not appear in the Stoner factors K dd (0, 0) 
and K c f(0,0). Moreover, we shall see below that the solutions to Kdd{0,0) = and 
K c f(0,0) = both lead to the same critical temperature. Therefore, the system's 
tendency for a certain Cooper pair symmetry only shows up in the intensity of A^O, 0) 
or A c j(0,0), which is controlled by the numerator of these functions. We also see that 
both antiferromagnetism and superconductivity are controlled by the same interaction 
parameter J, which, in turn, depends on hybridization only. 

In Figure [T] we show a plot of the Neel and superconducting temperatures as 
functions of the total electronic density n. It is seen that antiferromagnetism can 
only occur in a very small region of electronic density. Furthermore, the increase of 
T/v when n — > 2 corresponds to an increase of the density of n/ electrons towards 
the Kondo limit (rif — 1). It is also clear that T/v is not a monotonically increasing 
function of J. Upon reducing V, a larger range of electronic densities can be reached 
where antiferromagnetic order can be found. From the inset of Figure Q we see that the 
superconducting temperature T c is very small (about T c ~ T N /50) close to the density 
region where the system exhibits antiferromagnetic order. 

The dependence of the Neel and superconducting temperatures on pressure has 
been measured in some heavy-fermion systems jU EOl EI]- I n those studies the Neel 
temperature is found to decrease as the applied pressure increases and superconducting 
order is found to develop in a limited range of applied pressures, when the Neel 
temperature is reduced below ~1K. Let us now see how the critical temperatures in our 
model vary with the model parameters which, in principle, should depend on externally 
applied pressure. Increasing pressure should, presumably, make both the hybridization 
V and the conduction band hopping t increase [221 E3]- In Figure El we present T N and 
T c versus V, taking the ration V/t constant. We see that above a certain value of V 
the magnetic order disappears but the superconducting order remains. We also find 
that close to the region where the magnetic order vanishes, T c is much smaller than the 
maximum value attained by T/v- This is in qualitative agreement with experimental 
data on some Cerium compounds (e.g. Celn 3 ), where the ratio T N /T C ~ 100. Other 
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Figure 1. Neel Temperature as functions of the electronic density n. The inset 
shows the superconducting critical temperature as function of n. The temperatures 
are normalized by half bandwidth, D — 6t and t = 1. 



examples are: CeCu 2 (Sii_ a: Ge :2: ) 2 , where T c as function of pressure displays a positive 
curvature; and CeRhIn 5 , where the T c curve is almost parallel to pressure axis. [2E] 
For CeCu 2 Ge 2 and CeCu 2 Si 2 the T c curve initally stays almost parallel to the pressure 
axis, but it shoots up above a certain pressure. Although T c keeps increasing as V 
increases, it never reaches values comparable with the maximum value of Tjv, even for 
unreasonable values of V as we can see in the right panel of Figure 121 We believe that 
a better treatment of the boson propagator will lead to a decrease of T c in agreement 
with the experiments. We remark that the above calculation of T c is only valid in the in 
a situation where the system is non-magnetic because we have not calculated Kdd(0, 0) 
or K c f(0, 0) in the magnetically ordered phase. Furthermore, when T c is small, the 
approximation employed for the boson propagator should be improved by including its 
low energy part. 

For comparison we also plot the temperature T# defined as the difference between 
the renormalized /—level energy e/ and the chemical potential [HI EH lllj . This can 
be very different from the lattice Kondo temperature [21] in the non-magnetic system. 
Nevertheless, the combined behavior of T/v and Tk represents the well-known Doniach 
diagram showing the interplay between RKKY and Kondo screening effect. For small 
values of V, Tk is exponentially small and the system shows antiferromagnetic order. 
On the other hand, as V increases the Kondo temperature grows, leading to Kondo 
compensation of the /—moments and to a decreasing Neel temperature. For even 
larger values of V complete disappearance of the magnetic order takes place and the 
system shows paramagnetic behavior (assuming there are enough conduction electrons 
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Figure 2. Left: Neel Temperature T/v and superconducting critical temperature T c 
as function of V, for a constant ratio of V/t = 1.2, electronic density n = 1.8, and 
eq = —0.251? . Right: T c over a very large (nonphysical) values of V. Note that 
T c <t^. Tn always, (cc) and (cf) means that T c has been computed using equations 1221 
and 1231 respectively. Both give the same results. 



to compensate all the /—local moments). We have also computed the superconducting 
critical temperature from both equations and and obtained the same T c , as 
can be seen in the right panel of Figure EJ Along the Tjv curve, rif decreases from 1 to 
0.8, as V increases, and nj ~ 0.85 when T/v is maximum. 

Both Figure ^and Figure El show similar behaviour near the point where T/v — > 0. 
In both cases T c starts to increase with a positive curvature. Although in many heavy- 
fermion systems T c presents a negative curvature, there are examples where a positive 
curvature have been observed, such as CePd2Si2 under chemical pressure (T/v = 10 K, 
T c = 0.2 K, therefore T N /T C = 50) 

Since our treatment does not take competition between magnetism and 
superconductivity into account, we cannot tell whether finite values of T c and Tn 
imply that both types of order will be present at low temperature. Nevertheless, 
we found in previous work (2H1, at the simplest mean field level, that magnetism 
and superconductivity may coexist in the system. It follows from the above remarks 
that the calculation of T c when T/v is finite requires both the introduction of a better 
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approximation for the boson propagator and extra electronic propagators describing the 
antiferromagnetic order in the system, as was done in the description of spin waves in 
the magnetically ordered Mott insulator.) 
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